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Introduction 

Inventory  maintenance  is  very  critical  in  maintaining  fleet  readiness. 

There  is  a  cost  in  maintaining  a  set  inventory  and  a  cost  in  not  being  able 
to  supply  spare  parts  for  existing  assets.  Therefore,  the  Naval  Supply  Depot, 
under  Task  Order  MIPR2-0007,  requested  NAVAIRDEVCEN  to  develop  an  interactive 
computer  program  that  determines  if  there  is  a  statistically  significant 
positive  or  negative  trend  in  inventory  over  a  number  of  years.  A  negative 
trend  would  of  course  require  a  reordering  of  additional  spare  parts  while  a 
positive  trend  would  indicate  a  reduction  in  spare  parts,  thereby  reducing 
the  costs  of  inventory  maintenance. 

The  program  developed  and  discussed  in  this  report  utilizes  nonparametric 
statistics  to  test  if  the  inventory  changes  are  random  against  a  monotonic 
trend.  When  the  set  of  yearly  inventories  are  ranked  and  tested  against 
the  inventory  years,  the  agreement  or  disagreement  with  a  positive  yearly 
trend  determines  if  the  inventory  has  increased  or  decreased  over  the  years, 
i.e.,  a  large  positive  (negative)  number  would  indicate  a  positive  (negative) 
trend  while  a  number  close  to  zero  would  indicate  no  inventory  change. 

Mathematical  Background 

The  classical  measure  of  association  between  two  variables  is  the  correlation 
coefficient.  The  importance  of  the  correlation  coefficient  arises  from  para-  . 
metric  linear  regression  theory.  However,  in  non-normal  populations,  dependence 
among  variables  can  manifest  itself  in  other  ways.  For  example,  the  Kendall 
Rank  Correlation  coefficient  can  be  used  to  determine  the  degree  of  relation¬ 
ship  between  populations  and  is  dependent  on  the  value  of  parameter  "S".  The 
computation  of  S  is  performed  by  counting  how  often  the  two  rankings  from 
different  variables  move  in  the  same  direction  and  how  often  they  move  in  the 
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opposite  direction.  This  computation  is  performed  for  all  possible  pairs 
in  the  number  sequences.  The  strength  of  the  relationship  is  obtained  by 
considering  S  relative  to  its  two  extremes,  forming  the  quantity  "t"  known 
as  the  Kendall  Rank  Correlation  Coefficient.  A  value  of  t  near  +1  implies 
close  agreement  between  the  rankings  while  a  value  of  t  near  -1  implies 
opposite  rankings.  A  value  of  t  in  the  neighborhood  of  zero  indicates  neither 
agreement  nor  disagreement  (reference  (a)). 

Since  many  of  the  values  derived  for  t  and  consequently  S  do  not  differ 
significantly  from  zero,  a  statistical  test  must  be  applied  to  determine  when  ~~ 
t  differs  significantly  from  zero.  This  test  is  known  as  the  test  of 
"independence"  and  is  defined  in  reference  (a) . 

A  statistic  z  can  be  developed  to  test  the  hypothesis  Hq  that  the  two 
variables:  years  and  the  ranking  for  the  inventory,  are  independently  distributed, 
i.e.,  the  test  statistic  S  does  not  significantly  differ  from  zero  implying  that 
neither  a  positive  nor  negative  trend  exists.  The  test  statistic  z  for  a  sample 
with  no  ties  and  sample  size  of  at  least  eight  samples  can  be  approximated  by 
a  normal  distribution  and  is  given  by: 

*  6  S 

z  “  — — 

2n(n-l)(2n+5) 

The  program  as  defined  in  the  subroutine  TAU1  will  compensate  the  test 
statistic  z  for  ties  and  requires  a  minimum  of  six  samples. 

Program  Description 

The  trend  analysis  program  developed  is  an  interactive  program  written  to 
guide  the  user  through  the  program.  Following  the  sign-on  procedure  given  in 
Appendix  A  the  user  loads  the  program  from  the  library  with  the  command: 


)LOAD  TZIMES 
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The  trend  analysis  program  can  be  activated  by  the  coomtand  TREND  and  "return". 

The  program  will  then  request  the  user  to  enter  the  FSC  #  shorn  below 

enter  frc  mm  per 
Ps 

.  34 

The  #34  was  entered  and  the  program  responded  with 

rrc  supper  *34 

The  computer  then  requested  the  inventory  to  be  entered  for  each  year  by 

EE  TER  INVENTORY  PEI' AND  FROt:  FIRST  YEAR  TO  FREER  NT 

r  s 

The  inventory  495  through  550  was  then  entered  into  the  program  with  the 
depression  of  the  return  key 

495  525  520  490  555  530  475  510  515  545  540  550 

The  program  then  requested  the  user  to  enter  the  corresponding  years  by  printing 

ENTER  INVENTORY  YEAFE 
P : 

The  years  entered  oust  be  in  increasing  order  and  correspond  to  the  inventory 
entered.  For  this  problem  the  following  years  were  entered 

1970  1971  1972  1973  1974  1975  1970  1977  197P  1979  15CG  1901 

The  program  then  analyze*  the  entered  data  and  prints;  the  Kendall  rank  correlation 
coefficient,  the  S  value  and  the  descriptive  level  of  the  test,  where  this  level 
describes  the  area  under  the  normal  curve  to  the  right  of  the  computed  z  value 
defined  above.  The  smaller  this  value  is,  the  lower  the  probability  of  saying 
there  is  a  positive  trend  or  that  there  is  no  trend  (Type  II  Error).  The 
probability  of  a  type  two  error  is  8.5Z. 
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ccppelaticf  coefficifft=o . 303030303 

TPCPP  IF  POFITIVE  r=  2  0 

msCFIPTIVE  LEVEL  OF  TFFT  =  0.0f,:il7P?°9C 

The  program  then  requests  the  user  to  respond  to  the  following  cue 


u'ti.r  2  ir  you  vi nr  to  optaif  life  of  plcflcfiof  otefevih:  eptlf 


A  response  of  lf  skips  this  part  of  the  program  while  a  response  of  2  results 
in  the  following: 

1.  The  program  fits  a  straight  line  to  the  inventory  as  a  function  of  the 
entered  years.  This  program  is  listed  and  described  in  Appendix  B. 

2.  The  coefficients  A  and  B  of  the  linear  equation  Y  “  X  +  BX  are 


computed  and  listed. 

3.  The  standard  error  of  the  regressed  line  is  computed  and  listed. 

4.  The  correlation  coefficient  is  computed  and  listed. 


5.  The  T  value  from  which  the  confidence  can  be  computed.  These  outputs 


are  shown  below. 


■*  LE  FT  PAIC  UT  LIVE  IF  YsA*PX  VI™V  A  AFP  P  — — *>3I>0  407^25 
FT AFP ARP  cprrp  of  EFTIHATE  *0.4265687474 
CIVPLL  rOpFELATIOi:  COP FFI CIFP T  *0.42656374  74 
-m  VAT, UP  *1.491427002 


2.97202797 


The  program  then  plots  the  regressed  line  (circles)  and  the  inventory  (plusses). 
The  scatter  of  the  inventory  date  relative  to  the  regressed  line  indicates 
graphically  the  goodness  of  fit. 
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475.00000+  + 

/♦ - - - - - - - - - ♦ - + - + 

1970. 0C00  1971.8333  1973.  F.G6  7  1975.5000  1977.3333  1979.1607  1981. OOCO 


If  the  user  wishes  to  project  this  linear  trend  into  the  future,  the 
cue  response  to  "Inventory  Prediction"  should  be  "2"  and  return. 


rrrrj.  2  ir  ivvevtory  pfmictioi:  ir  fitiffp  crnzniin:  evtfp  1 

r* 

7 

The  program  then  requests  the  years  prediction  to  be  entered.  In  response 
to  the  entry  the  program  computes  and  lists  the  inventory  prediction  for  the 
entered  years. 


El' TER  YEARS  OF  rREVICTJCll 

p* 

1982  1983  1984 

IFVEF70FY  =540.1515152  543.1235431  540.0955711 
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The  program  can  then  be  cycled  or  ended  by  the  response  of  2  or  1  to  the 
next  inquiry. 

rmr  2  if  fcpl  fata  if  to  he  mat.itlf  cti  i.ih:  t.  eetep  1 
[•: 

o 

4. 

Summarv/Conclusion 

The  program  described  in  this  report  calculates  the  Kendall  Rank  Correlation 
coefficient  which  is  combined  with  a  test  for  independence  against  a  monotonies lly 
increasing  function.  This  problem  will  be  sensitive  to  changing  trends  when 
several  years  are  tested,  but  will  tend  to  be  insensitive  to  short  term 
inventory  changes. 
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APPENDIX  A 
Main  Program  -  TREND 


The  function  "TREND"  is  the  main  program  for  the  Trend  Analysis  program 
and  is  listed  below. 

7 TREND  CC]7 
7 TREND 

[1]  SP«-6  pO 

[ 2 ]  RE  TURK  1 : » ENTER  FSC  NUMBER ' 

[3]  ' FSC  NUMBER  =  •  ;[] 

[4]  • ENTER  INVENTORY  DEMAND  FROM  FIRST  YEAR  TO  PRESENT ' 

[  5 ]  SET+O 

[  6 ]  Y-SET 

[7]  'ENTER  INVENTORY  YEARS' 

[  8  ]  */?•*-□ 

[9]  -*•(  (pYR)*(pY)  )/ERROR 

[10]  ORDER  Y 

[11]  TAU1 

[12]  'CORRELATION  COEFFICIENT-' ; ( 2*S) * (P* (P-1 ) ) 

[13]  +(S<0)/PROB 

[14]  'TREND  IS  POSITIVE  S=';S 

[15]  'DESCRIPTIVE  LEVEL  OF  TEST  = ' ; 1 - ( NCUM  Z) 

[16]  RETURN2:' ENTER  2  IF  YOU  WISH  TO  OBTAIN  LINE  OF  REGRESSION  OTHERWISE  ENTEl 

[17]  URNo-U 

[18]  ■+•  ( URN- 1 )  /  RETURN3 

[19]  T*-YR  SR  SET 

[20]  £IP/?rG«-(r[4;l]xy/?)+r[3  jl] 

[21]  'THE  STRAIGHT  LINE  IS  Y~A*BX  WITH  A  AND  B  = ' ; T[ 3 ; 1 ] ,TL 4 ; 1 ] 

[22]  • STANDARD  ERROR  OF  ESTIMATE  =*;2’[5;1] 

[23]  'SIMPLE  CORRELATION  COEFFICIENT  =  *  ;  7*[  5  ;  2  ] 

[24]  ' T  VALUE  *  *  s T[ 4 ; 3 ] 

[25]  MAXLIN+l/l  l/LINREG  J /SET) 

[26]  MINLIN+l / l l/LINREG ,1/SET) 

[27]  SPl  1  ]«-30 

[28]  SPl2]*-MINLIN 

[29]  SPlSl*-(MAXLIN-MINLIN)*S 

[30]  SPl  4  ]«-6  0 

[31]  SP[5]«-L/yP 

[32]  SPl  6]«-(r/y/?-L  /YR)  *6 

[33]  SP  PLOT  LINREG  AND  SET  VS  YR 

[34]  'ENTER  2  IF  INVENTORY  PREDICTION  IS  DESIRED  OTHERWISE  ENTER  1' 

[35]  URN 2*-D 

[36]  -*(UFN2S1  )/  RETURNS 

[37]  ' ENTER  YEARS  OF  PREDICTION ' 

[38]  PREY*-  0 

[39]  'INVENTORY  ='  ;  (  T[  4  ;  1  ]  xPPPr ) +r[  3  ;  1  ] 

[40]  RETURNS -.'ENTER  2  IF  MORE  DATA  IS  TO  BE  ANALYZED  OTHERWISE  ENTER  1' 

[41]  URN1+V 

[42]  -*(URN  l=l)/0 

[43]  ^RETURN  1 

*44]  ERROR-.' INCONSISTENCY  BETWEEN  DATA  BASE  AND  YEARS' 

■45]  *-RETURN  1 

[46]  PROBx  'TREND  IS  NEGATIVE  S=  ' iS 

47]  'DESCRIPTIVE  LEVEL  OF  TEST  * ’ ; 1- (NCUM \ Z ) 

48]  +RETURN2 
7 
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FUNCTION 

Shapes  a  row  vector  SP 

Requests  user  to  enter  FSC  # 

Prints  entered  FSC  # 

Requests  inventory 

Assigns  the  variable  SET  equal  to 
the  inventory 

Assigns  the  variable  Y  equal  to  SET 

Requests  inventory  years 

Assigns  the  variable  YR  to  inventory 
years 

Checks  that  each  year  has  an  inventory  - 
if  there  is  an  inconsistency,  a  skip  to 
statement  M  occurs 

The  function  ORDER  is  called  and  the 
rank  for  the  inventory  is  determined 
as  a  function  of  the  years 

The  function  TAUl  is  then  called  and  the 
variable  S  is  determined 

Kendall  Rank  correlation  coefficient  is 
computed 

"SIGN"  of  S  is  determined  and  if  negative 
skipi  to  statement  [463 

Positive  value  of  S  is  noted  and  value 
is  printed 

Descriptive  level  of  test  -  (area  under 
normal  curve  to  the  right  of  "Z") 

Open  the  option  of  user  to  obtain 
coefficients  for  linear  equation  and 
obtain  for  plot 

URN  -  is  the  response  variable  of 
statement  Cl6l 

Eest  of  statement  with  skip  to  statement 
40 J  if  one  is  entered 


A2 
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STATEMENT  # 


[27  -  [32! 


FUNCTION 

Call  for  regression  SR  Chat  computes 
coefficients  and  fits  statistics 

LINREG  is  the  linear  equation  for  the  . 
line  with  coefficients  computed  in  19 


Coefficients  A  and  B  are  printed  for 
curve  Y  *  A  +  BX 

The  standard  error  of  the  estimate  is 
printed 

The  correlation  coefficient  is  printed 

The  T  value  is  printed 

Maximum  value  of  the  regressed  curve  is 
determined 

Minimum  value  is  determined  for  regressed 


Scaling  is  performed  for  the  plot  function 

The  plotting  function  is  called  and  the 
data  vs  Y  are  plotted 

Prediction  is  selected  by  user  or  program 
skips  to  ^0 

URN  2  is  the  variable  name  for  entry  into 
prediction 

The  variable  URN  2  is  tested  for  1  or  2 

Request  for  years  to  which  prediction 
is  desired 

PREY  -  variable  nasie  for  the  prediction 
years 

Inventory  is  computed  and  outputted 

Enter  new  set  of  data  or  end  program 

Variable  name  URN  1  for  response  to 
item  ^(j 

Test  of  variable  URN*  program  goes 

to  logical  zero  which  is  end  of  program 

Go  to  statement  #  ! 2] 


A3 
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STATEMENT  # 

H 

H 

H 

H 


FUNCTION 

Go  to  statement  when  there  is  inconsistency 
between  YR  and  SET  variable 

Go  to  statement  #  ^2j 

Go  to  statement  when  S  is  negative 

Descriptive  level  of  test  for  negative  Z 

Go  to  statement  #  [2^ 
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APPENDIX  B 
Subroutine  -  ORDER 


This  function  ranks  the  data  entered  into  the  program  and  counts  the 
number  of  ties.  The  variable  to  be  ranked  is  X. 


VOPTEF  CT]V 
VOPPEF  X ; V 
1-1-3 - K-l - 

[2]  TIL+C 

[3]  y<-ox(ipA') 


[  4  3 

rr„rnf,'  y=[/v)  /V 

is] 

■*(  (pV)>2) /ADJUST 

[G] 

r<-x  1  v 

-C-7-] —  . 

-lUH' - 

[83 

A'[P]«-9  99  ■  •  •  r: 

[93 

+ 1 

ri-C  3 

-*  (  y  >  p  Y  )  /C 

ill] 

■*.petupi: 

[12] 

ADJUST:  L+ 1 

:m- 

— TESTIS*  1 - , 

[14] 

LETUP!' 1  :  F*-X  1 1'[  1  ] 

[19] 

•:i«-l*  (  pT/) 

C  1C  3 

iC  ?  3  i 

[17] 

y[  r  ]4-°  99 

[18] 

l+L*  1 

[19] 

+(L>pV)/ADJUSTl 

[20] 

-*RETURVl 

[21] 

AEJVST1  :/>/>  (  pT'> 

r  22  1 

>£>y  )/c 

[23] 

-PETURD 

STATEMENT  # 


FUNCTION 


Initialization  of  index  K 

Initialization  of  the  variable  TIE 

Initialization  of  the  row  vector  Y 
with  length  of  the  input  variable  X 

V  is  a  row  vector  containing  the 
minimum  elements  of  the  vector  X 


*  n 
5; 


H 


V  contains  two  values  i.e.,  there 
.a  ti^  the  program  jumps  to  statement 


Bl 
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STATEMENT  # 


FUNCTION 


i6. 


P  is  index  variable  for  the  element  V 
in  the  vector  X 


The  element  P  of  vector  Y  is  set  equal 
to  K 


The  element  P  of  vector  X  is  set  equal 
to  a  large  number 


The  index  K  is  incremented 

End  program  if  the  index  K  is  greater  than 
the  #  of  elements  in  X 


H- 

-  h 

H 


(l3 


Go  to  statement  ^4j 

The  subroutine  for  ties  within  the  matrix 
elements 

Set  variable  L  equal  to  one 
TIE  is  incremented 


Set  the  index  P  to  the  first  element  of  V 


K1  is  the  fraction  to  be  added  to  K  for 
the  rank  if  more  than  one  element  are 
ranked  the  same 


The  rank  of  yLp}  is  set 


L  5 

H 

H 


The  variable  X[Pj  is  set  to  a  large  # 
Increment  the  variable  L 

Go  to  statement  .  [2  J  if  all  of  the  elements 
that  are  equal  are  set  to  the  same  rank 


H 

H 

N 


Go  to  statement 


Set  the  rank  #  K  to  the  next  value  which 
equals  the  last  rank  plus  ties 


End  program  if  all  elements  in  X  are 
complete 

Go  to  statement 
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APPENDIX  C 
Subroutine  -  TAU1 


The  variable  S,  calculated  by  the  subroutine*  determines  the  number  of 
times  the  ranking  moves  in  the  same  or  opposite  direction  from  the  monotomic 
function. 


VTAV1  [n]V 
VTAL1 
[ 1 ]  TU-2 

C  2  ]  Z>0 

[  3  ]  F+-0 

[4]  A«-pY 

[  5  ]  A>0 

tc]  r:  rix-c 

[7]  A?Cl«-0 

[ 8 ]  BCTURBtK+Z* 1 

[  9  ]  F+r.  +  iUt-K) 

tie]  n+pp 

[11]  rc*-+/(  (RpYlKi  )<y[p] ) 

[12]  ATW  ((Tpy[^])>J[P]) 

[13]  FP1+RD1+KD 

[14]  A  Cl  «-ArCl  ♦A?C 

[15]  s+s+trc-ru) 

[16]  ■+(/f=  A')  /COH 

[17]  •*  PH  TUP  I' 

[15]  CCA':  VAFSl<-!t*  (  A  - 1  )  *  (  (  2  *  A' )  +  5  ) 

[19]  -*■(  ?in<  i  )/coi;fipepce 

[20]  AtJVSTiL+L+{?U*{?U-l)*{  (2*TV)*S) ) 

[21]  ■♦(  OU-TIE)  / COEFItEVCE 

[22]  TV+7U+ 1 

[23]  -+ADJVST 

[24]  COF  FI  DEFOE :  (  F4.»?S1 -£ )  *  1  8 

[25]  SICS-VAPS*  0.5 

[26]  C-StCICS 

[27]  *0 
V 


STATEMENT  # 


»1 


l4i 


FUNCTION 

Initialize  variable  TU 

Initialize  variable  L 

Initialize  variable  S 

Set  the  variable  N  equal  to  number  of 
element  in  the  row  vector  Y 


Cl 
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* 


« 


I 


STATEMENT  # 

V 

'[*) 

M 

w 

w 

H 

H 

H 

H 

H 

H 

H 

H 

4 

H 

H 

H 

H 

r  m 


FUNCTION 


Initialize  variable  K 

Initialize  variable  ND1 

Initialize  variable  NCI 

Increment  variable  K 

Set  P  as  the  monotomic  function  with 
the  same  number  of  elements  as  X 

R  is  the  number  of  elements  in  P 

Compute  the  number  of  elements  that 

increase  monotonically  in  Y  i.e.,  are 
concordant 

Compute  the  number  of  elements  that 
decrease  monotonically  in  Y  i.e.,  are 
discordant 

Counter  for  the  number  of  discordant  values 
Counter  for  the  number  of  concordant  values 


S  is  the  total  difference  between  concor¬ 
dant  and  discordant  values 


If  the  index  K  is  equal  to  the  nunjp^r 
of  elements  in  Y  go  to  statement 


Go  to  statement 


M 


impei 

ft 


Compute  the  variance  of  S 

If  no  ties  in  the  ranks  go  to  statement 

Compute  a  correction  for  the  variance  of 
S  for  ties 

If  the  variable  TU  equals  the  number  of 
ties  go  to  statement  ^ 

Increment  the  variable  TU 


Go  to  statement 


Compute  the  variance  of  S 


Compute  the  standard  deviation  of  S 
Compute  the  statistic  Z  with  N(0,1) 
End  program 


C2 


J 
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APPENDIX  D 
Subroutine  -  NCUM  2 


This  subroutine  is  an  approximation  for  the  inverse  normal  distribution. 


7 ;;,r" 


[1]  ?rp+-  ( 4 .  3  <■'  6  3  BT~  5  *  Z  *6  )  ♦  (  C  .  00  C  2  7  B  5  £  7  2  *Z  *5  )  +  (  C  .  C  0  0 1  5  2  0 1 4  3  *  Z  *  4  ) 

[2] *C"  rrof*^i-i*r/“p*is‘C422£2ri22,<2*2)  +  (^ .n70523r7?2?799e*z)*l 


D1 
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APPENDIX  E 
Subroutine  -  SR 


Subroutine  SR  ia  pert  of  the  acientific  aubroutine  package  in  the  APL 
library.  The  following  documentation  ia  available  for  thia  program.  A 
description:  SRHOW  and  a  program  liating  SR. 


SRHOV 

SIMPLE  REGRESSION 
T+X  SR  I 
ENETREDi 04/16/70 

X  AND  Y  ARE  VECTORS  GIVING  THE  ( SAME  NUMBER  OF) 
'OBSERVATIONS  ON  AN  INDEPENDENT  VARIABLE  X  AND  A 
DEPENDENT  VARIABLE  I.  T  IS  A  MATRIX  WITH  5  ROWS 
AND  3  COLUMNS  CONTAINING  THE  RESULT  OF  FITTING 
THE  STRAIGHT  LINE  Y=A+B*X  BY  THE  METHOD  OF  LEAST 
SQUARES  IN  THE  FOLLOWING  FORMAT : 

ROWli  MEAN  OF  X,  ST  DEV  OF  X .  0 
ROW2i  MEAN  OF  Y .  ST  DEV  OF  7.  0 
ROWS-.  A,  0.  0 

ROWki  B.  ST  ERROR  OF  Bt  T- VALUE 

ROWS :  ST  ERROR  OF  ESTIMATE .  R*SIMPLE  CORR  COEFF .  R* 2 


VSR  ED3 V 
VT+X  SR  J 

[1]  SX<-(  U++/{X»MX+{+/X)*n)*2  )♦(»«-(  px)  )-l  )*0.5 

[2]  SY+(  (B'»+/(Y-MY+(+/Y)*N)*2)*N-l  )*0.5 

[3]  BO+MY-MX*Bl+[+/{X-MX)*(Y-MY) )«A 

[4]  SE+(  (B*l-RS<J+(R+Bl*SX*SY)*2)*N-2)  +  0.S 

[5]  TV+B1*SB1+(SY*SX)*( (N-2 )*(1-RSQ) )*0.5 

C  6 ]  JV(5  3)pMXtSX,0,MYtSYt0tB0,0  0 ,B1 ,SB1 .TV ,SE ,R ,RSQ 
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APPENDIX  F 
Subroutine  -  PLOT 


The  plotting  routine  i«  also  part  of  the  aceintific  subroutine  in  the 
APL  library.  The  Plothov  function  describes  its  use  and  capability  and  the 
prograa  is  listed. 


PLOTHOV 

CRAPHING  FUNCTION :  C+S  PLOT  D 

DEPENDENT  VARIABLES  ARE  PLOTTED  AGAINST  AN  INDEPEN¬ 
DENT  VARIABLE.  IF  THE  RIGHT  ARGUMENT  D  IS  A  VECTOR  OP  A 
l -ROW  MATRIX  THE  ELEMENTS  OF  D  ARE  CONSIDERED  TO  BE  VALUES 
OF  A  DEPENDENT  VARIABLE.  AND  THE  VALUES  ARE  PLOTTED  AGAINST 
THIER  INDICES.  IF  D  IS  Z  MATRIX  WITH  >  1  ROW,  THE  FIRST  ROW 
SHOULD  CONTAIN  THE  VALUES  OF  THE  INDEPENDENT  VARIABLE  AND 
THE  REMAINING  ROWS  SHOULD  CL  TAIN  CORRESPONDING  VALUES  FOR 
THE  DEPENDENT  VARIABLES. 

THE  LEFT  ARGUMENT  S  INDICATES  SCALING.  IF  OMITTED  THE 
GRAPH  WILL  BE  SCALED  AUTOMATIC ALLI  TO  BE  40  TYPEWRITER  LINES 
HIGH  BY  70  TYPEWRITER  SPACES  WIDE.  IF  YOU  WISH  TO  SPECIFY 
THE  SIZE  ONLY  AND  ALLOW  THE  PROGRAM  TO  SCALE  AUTOMATICALLY. 

THE  VALUES  OF  S  SHOULD  BE  AS  FOLLOWS: 

SCI]  =  THE  DESIRED  HEIGHT  OF  THE  GRAPH  IN  LINES  (AT  LEAST  5) 
S[ 2 ]  =  THE  DESIRED  WIDTH  OF  TEE  GRAPH  IN  SPACES  (AT  LEAST  10) 

IF  YOU  PREFER  TO  SPECIFY  THE  SCALING  OF  THE  GRAPH 
YOURSELF.  YOU  MUST  GIVE  IN  ADDITION  THE  VALUES  FOR 
THE  SCALE  ORIGIN  AND  THE  DESIRED  DIFFERENCES  BETWEEN 
SCALE  MARKS.  IN  THIS  CASE  S  SHOULD  BE  ASSIGNED 
VALUES  AS  FOLLOWS: 

SCI]  =  THE  DESIRED  HEIGHT  IN  LINES  (AT  LEAST  5) 

SC  2]  *  THE  SCALE  ORIGIN  FOR  THE  VERTICAL  AXIS 
SC  3]  =  THE  DIFFERENCE  BETWEEN  SCALE  MARKS  ON  THE 

VERTICAL  AXIS  (SCALE  MARKS  ARE  EVERY  5  LINES ) 

SC 4]  =  THE  DESIRED  WIDTH  IN  SPACES  (AT  LEAST  10) 

SC  5]  =  THE  SCALE  ORIGIN  FOR  THE  HORIZONTAL  AXIS 
SC 6]  =  THE  DIFFERENCE  BETWEEN  SCALE  MARKS  ON  THE 

HORIZONTAL  AXIS  (SCALE  MARKS  ARE  EVERY  10  SPACES ) 


IF  SEVERAL  DEPENDENT  VARIABLES  ARE  PLOTTED.  THEY 
ARE  SHOWN  ON  THE  GRAPH  USING  SUCCESSIVE  SYMBOLS  FROM 
THE  VARIABLE  PS.  IF  PS  DOES  NOT  EXIST.  IT  WILL  BE  CREATED 
WITH  THE  VALUE"' o*D*lTA7«'.  IF  THE  PLOT  SYMBOLS  IN  PS  ARE 
EXHAUSTED.  THE  FIRST  ARE  USED  AGAIN.  THE  VARIABLE  PjrMAY 
BE  REASSIGNED  TO  CAUSE  DIFFERENT  PLOT  SYMBOLS  TO  BE  USED. 

IF  TWO  OR  MORE  VARIABLES  HAVE  A  VALUE  IN  THE  SAME  PRINT  POSITION . 
THEIR  SYMBOLS  WILL  BE  OVERSTRUCK  AT  THAT  POSITION. 


1 


NADC-82246-20 


IP  THE  VARIABLES  DX  AND  DJ  ARE  DEFINED  TO  BE 
CEARACTER  ARRAYS  CONTAlTTNG  DESCRIPTIONS  OF  THE 
INDEPENDENT  AND  DEPENDENT  VARIABLES  RESPECTIVELY .  THE 
DESCRIPTIONS  WILL  BE  PRINTED  ALONG  THE  APPROPRIATE  AXES. 
IF  THESE  ARE  EMPTY  VECTORS  OR  DO  NOT  EXIST 
NO  DESCRIPTION  WILL  BE  PRINTED. 

THE  RESULT  G  IS  A  CHARACTER  VECTOR  WITH  EMBEDDED 
CARRIAGE  RETURNS.  IT  APPEARS  TO  BE  A  MATRIX  WITH  SHAPE : 
SIZE* 4  13 


WHERE  SIZE  IS  TEE  2 -ELEMENT  VECTOR  SPECIFYING  HEIGHT .WIDTH 
OF  THE  PLOT  IN  LINES .SPACES  (40  70  IP  S  NOT  SPECIFIED). 

A  FUNCTION  ' VTOM '  IS  PROVIDED  TO  CONVERT  SUCH  A  VECTOR 
RESULT  TO  A  MATRIX  WITH  THE  SAME  APPEARANCE . 

TEE  FUNCTIONS  ' VS '  AND  'AND'  ARE  PROVIDED  TO  HELP 
FORM  THE  RIGHT  ARGUMENT  TO  PLOT .  THESE  FUNCTIONS  ALLOW  VECTOR 
OR  MATRIX  ARGUMENTS  AND  FORM  A  MATRIX  RESULT  BY  JOINING  THE  ROWS 
IN  THE  APPROPRIATE  ORDER.  FOR  EXAMPLE. 

A  AND  B  AND  C  VS  D 

WOULD  JOIN  THE  VECTORS  A.  B.  C.  AND  D  TOGETHER  TO  FORM  A  MATRIX 
HAVING  ROWS  TAKEN  FROM  D,  A.  B,  AND  C.  IN  THAT  ORDER. 

THERE  IS  SUFFICIENT  WORKING  AREA  TO  EXECUTE  IF 
DWA  Cl]  i  (  0 . 4  **  /SJZD  r  (3  **/pP) 

ANOTHER  1100  WORDS  CAN  BE  RELEASED  BY  USING 
) ERASE  . GRPDOC 


EXAMPLES-. 


10  20  PLOT  ilO 
10* 


o 


o 


o 


o 


10 


F2 
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10  20  PLOT  (O.ilO)  AND  (11-tll)  AND  (llp5)  VS  C 
10  +  *  o 
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